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For every positive integer n let S(n) be the minimal positive integer m such 
that n | m!. For any positive number z > 1 let 


A(z) = = > S(n) - 


nz 
be the average value of S on the interval {1, z]. In [6], the authors show that 
A(z) <cr+c2 (2) 


where c, can be made rather small provided that z is enough large (for example, 
one can take cy = .215 and cg = 45.15 provided that z > 1470). It is interesting 
to mention that by using the method outlined in [6], one gets smaller and smaller 
values of c; for which (2) holds provided that z is large, but at the cost of increasing 
co! In the same paper, the authors ask whether it can be shown that 


2z 
A 
@)<m (3) 
and conjecture that, in fact, the stronger version 
Als} (4) 


logz 


might hold (the authors of [6] claim that (4) has been tested by Ibstedt in the 
range z < 5- 10° in [4]. Although I have read [4] carefully, I found no trace of the 
aforementioned computation!). 


In this note, we show that is indeed the correct order of magnitude of 


A(z). 
For any positive real number z let a(x) be the number of prime numbers less 
then or equal to z, 


x 
logz 


B(z)=zA(z)= QD) S(n), (5) 
E(z) = 2.5 loglog(z) + 6.2 + ~ (6) 
We have the following result: 
Theorem. 
5(a(2) — (/z)) < A(z) < a(x) + E(x) for all x > 3. (7) 


Inequalities (7), combined with the prime number theorem, assert that 
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log z logz 


which says that ae is indeed the right order of magnitude of A(z). The natural 
conjecture is that, in fact, 


A(z) = — +0( cd ). (8) 


~ jogz log’ x 
Since 


zr 1 z 3 
=. SS =< = > 
log z Ct Sale) S Geet cies? fore 2 08 


it follows, by our theorem, that the upper bound on A(z) is indeed of the type (8). 
Unfortunately, we have not succeeded in finding a lower bound of the type (8) for 
A(z). 


The Proof 


We begin with the following observation: 


Lemma. 
Suppose that n = p%'...pp* is the decomposition of n in prime factors (we 
assume that the p;’s are distinct but not necessarily ordered). Then: 


1. 
S(n) < maxy_,(ai7%). (9) 


2. Assume that a,p, = max!_,(aip;). If a; < pi, then S(n) = aipi. 
3. 
S(n)>ai(pi-1)  foralli=1, ..., &. (10) 


Proof. 

For every prime number p and positive integer k let e,(k) be the exponent at 
which p appears in k!. 

1. Let m > maxi, (aipi). Then 


€p,(m) = >, |= > [=| >a; fori=1, ..., k. 


oy : 
s>1 P; Pi 


This obviously implies n | m!, hence m > S(n). 

2. Assume that a; < pr. In this case, S(n) > a p;. By 1 above, it follows that 
in fact S(n) = ayp1. 

3. Let m = S(n). The asserted inequality follows from 


foo] 
m 1 m 
ai Se,(m) =o [=| <m > o= a 
t 


s>i * sp? 
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The Proof of the Theorem. 


In what follows p denotes a prime. We assume zr > 1. The idea behind the 
proof is to find good bounds on the expression 


B(z)- BiVz)= > S(n). (11) 


Jrinks 
Consider the following three subsets of the interval J = (./z, z]: 


C, = {n € I | S(n) ‘is not a prime}, 

Cr = {n€T| Sin) =p < va}, 

C3={ne€I| S(n)=p> Vz}. 
Certainly, the three subsets above are, in general, not disjoint but their union covers 
I. Let 


D(z) = > S(n) fori =1, 2, 3. 
ne€C, 


Clearly, 
max(D;(z) | i= 1, 2, 3) < B(z) - B(Vz) < Di(z) + Da(z) + Da(z). (12) 


We now bound each D; separately. 


The bound for D;. 


Assume that m € C,. By the Lemma, it follows that S(m) < ap for some 
p* || mand a@ > 1. First of all, notice that S(m) < a./m. Indeed, this follows from 
the fact that 
S(m) <ap<apt?<a/m  fora>2. 


In particular, from the above inequality it follows that p< /m < Vz. Write now 
m = p*k. Since m < 2, it follows that k < 2/p%. These considerations show that 


= z — a 2p—1 

Diz)< Yo or eat iat? gay 8) 
pS Vz o22 pS Vzra22 PSVE 

In the above formula (13), we used the fact that 


ae t=F(-)-15 (4) 1575 for |z} <1 


with z = 1/p. Since 


2p—1 5 
< — for p > 3, 
(p—1)? ~ 4p : 
it follows that 
1 1 
D(z) < 2(3—2 2 y -) = z(2.375 + 1.25 > =). (14) 
8 4 P P 
pSv= PS Vt 
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From a formula from [5], we know that 
1 
So = <loglogy+1.27 for ally>1.. 
psy 


Hence, inequality (14) implies 
Dy(z) < 2(2.375 + 1.25 (loglog Vz + 1.27)) <2(3.1+ 1.25 loglogz). (15) 


The bound for D, 
Assume that S(m) = p. Then m = py where p does not divide y. Since 
m > 1/2, it follows that 


r 
Pp P 
Since p < ,/z, it follows that at least one integer in the above interval is a multiple 
of p; hence, cannot be an acceptable value for y. This shows that there are at most 


ree 
Pp 2 P 
possible values for y. Hence, 
z- Vz 
Do(z) < D> p-( ) < (2- vz)x(vz). (16) 
P 
PSVz 
Bounds for D3 
Assume S(m) = p for some p > \/z. Then, m = py for some y < z/p. Hence, 
z 
D3(z)= D> p- RE (17) 


Vi<pse 


Notice that, unlike in the previous cases, (17) is in fact an equality. Since z > [z} > 
.5z for all real numbers z > 1, it follows, from formula (17), that 


52(n(2) — *(V2)) < Ds(z) < 2(x(2) — (V2). (18) 


Denote now by 
F(z) = 3.1+ 1.25 log log(z) 


From inequalities (12), (15), (16) and (17), it follows that 
52(a(z) — 1(/z)) < D3(z) < B(z) — B(/z) < Di(z) + Do(z) + Da(z) < 


2F (2) + (2 — V2)x(/z) + 2(x(z) — a(Vz)) = 24 (2) — Ven(V2) + 2F (2). (19) 


The left inequality (7) is now obvious since 


B(z) > B(/z) + .52(a(z) — a(/z)) > 1+ 52(a(z) — 2(V/z). 
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For the right inequality (7), let G(z) = x(x). Formula (19) can be rewritten as 
B(z) — B(/z) < G(z) - G(Vz) + 2F (2). (20) 
Applying inequality (20) with z replaced by /z, z'/4, ..., 21/?" until 2/2" < 2 
and summing up all these inequalities one gets 
bd . . 
B(z) — B(1) < G(z) + $02? F(2¥//). (21) 
s=0 
The function F(z) is obviously increasing. Hence, 
s Z 
B(z) <1+G(z)+ F(z) > 2”. (22) 
i=0 


To finish the argument, we show that 
z> Soa. (23) 
t=1 


Proceed by induction on s. If s = Q, there is nothing to prove. If s = 1, this just 
says that z > ./z which is obvious. Finally, if s > 2, it follows that z > 4. In 
particular, c > 2,/z or z — /z > \/z. Rewriting inequality (23) as 


s 


bm es 


which is precisely inequality (23) for ,/z. This completes the induction step. Via 
inequality (23), inequality (22) implies 


Biz) <1+2n(z) 4+ 22 F(z) = 14+ ea(z) + 22( 3.14 1.25 log logz 24) 
} 


or 
A(z) < #(z) + + +6.2+ 2.5 loglogz = x(x) + E(z). 


Applications 


From the theorem, it follows easily that for every ¢ > 0 there exists zg such 


that 
x 


A(z) < (14+ eee 


(25) 


In practice, finding a lower bound on Zo for a given €, one simply uses the theorem 
and the estimate 


xe 3 
n(z) < fogz (! + ae) for z > 1. (26) 
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(see [5]). By (7) and (26), it now follows that (25) is satisfied provided that 


z 1 3 
log z € Cre. r it E(2)) : 
For example, when ¢ = 1, one gets 
zr 
A(z) < aaae for z > 64, (27) 
for ¢ = .5, one gets 
z 
o-— > 
A(z) <1 Pines for z > 254 (28) 
and for « = 0.1 one gets 
A(z) <1.1—— for z > 3298109. (29) 
logz 


Of course, inequalites (27)-(29) may hold even below the smallest values shown 
above but this needs to be checked computationally. 
In the same spirit, by using the theorem and the estimation 


z 1 
—— a > 
n(z)> eae (1 + ard for z > 59 


(see (5]) one can compute, for any given ¢, an initial value zo such that 
A(z) > (.5-€) = forzr >z 
; logz e 


For example, when ¢ = 1/6 one gets 


——  forz > 59. (30) 


A(z) > pee 


ole 


Inequality (30) above is better than the inequality appearing on page 62 in (2] which 
asserts that for every a > 0 there exists rg such that 


A(z) > 2°/* for z > zo (31) 
because the right side of (31) is bounded and the right side of (30) isn’t! 
A diophantine equation 


In this section we present an application to a diophantine equation. The ap- 
plication is not of the theorem per se, but rather of the counting method used to 
prove the theorem. 

Since S is defined in terms of factorials, it seems natural to ask how often the 
product S(1) -.$(2)-...-.$(n) happens to be a factorial. 

Proposition. 


The only solutions of 


S(1) - $(2)-...-S(n) =m! (32) 
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are given byn=m€ {l, 2, ..., 5}. 

Proof. 

We show that the given equation has no solutions for n > 50. Assume that this 
is not so. Let P be the largest prime number smaller than n. By Tchebysheff’s 
theorem, we know that P > n/2. Since S(P) = P, it follows that P | m!. In 
particular, P < m. Hence, m > n/2. 

We now compute an upper bound for the order of 2 in S(1) - S(2)-...- S(n). 
Fix some § > 1 and assume that k is such that 2° |j S(k). Since 


S(k) = max(S(p*) | p* || £), 


it follows that 2° || S(p%) for some p® || k. 
We distinguish two situations: 


Case 1. 

p is odd. In this case, 2p | S(p*). If 8 = 1, then a = 2. If @ = 2, then a = 4. 
For 6 > 3, one can easily check that a > 29 —8+1 (indeed, if a < 28 —@, then one 
can check that p® | (2°p— 1)! which contradicts the definition of S). In particular, 
p?’-+1 | &. Since 27-1 > 2 +1 for z > 3, it follows that a > 28-142. Since 
k <n, the above arguments show that there are at most 

n 

28 


for 8 = 1, 2 


as) 


and 
pete for B > 3 

integers k in the interval [1, n] for which p | k, S(k) = S(p%), where a@ is such that 
p* || & and 29 || S(k). 

Case 2. 

p=2. If 8 = 1, then k = 2. If 8 = 2, then k = 4. Assume now that 6 > 3. 
By an argument similar to the one employed at Case 1, one gets in this case that 
a > 28 — B. Since 2% || k, it follows that 92°-8 | k. Since k < n, it follows that 


there are at most . 


Q3°=B 
such k’s. 


From the above anaysis, it follows that the order at which 2 divides S(1)- S(2)- 
...- S(n) is at most 


1 2 
a<itn (S444 oem) t+" pov (38) 
yeaa 28 b23 


(the number 3 in the above formula counts the contributions of S(2) = 2 and 
S(4) = 4). We now bound each one of the two sums above. 
For fixed p, one has 
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Hence, 


1 
» (4+5 ++ yp pe fas) < <L PoP oF 5 < 245 (40) 


We now bound the second sum: 


B 3. 4 5 B 
De g2°-6 «95 + 53 + p37 + 96 + » 92+4(8-2) 
p>3 
3 1 y+2\ 3.1715 , 31 
tad ror) = get gag aap) <0 a) 
From inequalities (38), (40) and (41), it follows that 
ep <3+.344n. (42) 


We now compute a lower bound for eg. Since e2 = e2(m!), it follows, from Lemme 
1 in [1] and from the fact that m > n/2, that 
log(m + Is n  log(n/2+ 1) 


ts gi eee a 
eames er ae log 2 (43) 


From inequalities (42) and (43), it follows that 


log(.5n + 1) 


. > .on— 
3+ .344n > .5n ieee 


which gives n < 50. One can now compute S(1) - S(2)-...- S(n) for all n < 50 
to conclude that the only instances when these products are factorials are n = 
Be 2 cod De 


We conclude suggesting the following problem: 


Problem. 
Find all positive integers n such that S(1), S(2), ..., S(n*) can be arranged in 
a latin square. 


The above problem appeared as Problem 24 in SNJ 9, (1994) but the range 
of solutions was restricted to {2, 3, 4, 5, 7, 8, 10}. The published solution was 
based on the simple observation that the sum of all entries in an n x n latin square 
has to be a multiple of n. By computing the sums B(z?) for z in the above range, 
one concluded that B(z?) # 0 (mod z) which meant that there is no solution for 
such z’ses. It is unlikely that this argument can be extended to cover the general 
case. One should notice that from our theorem, it follows that if a solution exists 
for some n > I, then the size of the common sums of all entries belonging to the 
same row (or column) is © nz(n?). 


Addendum 


After this paper was written, it was pointed out to us by an annonymous referee 
that Finch [3] proved recently a much stronger statement, namely that 


lim est) . A(z) = ~ = 0.82246703... (44) 


r—coO 
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Finch’s result is better than our result which only shows that the limsup of the 
expression log(z)A(z)/z when z goes to infinity is in the interval [0.5, 1]. 
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